METHOD FOR LOAD BALANCING A LINE OF 
PARALLEL PROCESSING ELEMENTS 



CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] The present application is related to U.S. Patent Application Serial No. _j 

entitled "Method for Using Extrema to Load Balance a Loop of Parallel Processing Elements" 

filed (DB00 1067-000, Micron no. 02-1297), U.S. Patent Application Serial 

No. entitled "Method for Load Balancing a Loop of Parallel Processing Elements" 

filed' (DB00 1068-000, Micron no. 02-1298), U.S. Patent Application Serial 

No. entitled "Method for Load Balancing an N-Dimensional Array of Parallel 

Processing Elements" filed (DB00 1062-000, Micron no. 02-1295), U.S. Patent 

Application Serial No. entitled "Method for Using Filtering to Load Balance a Loop 

of Parallel Processing Elements" filed (DB00 106 1-000, Micron no. 02-1296), U.S. 

Patent Application Serial No. entitled "Method for Rounding Values for a Plurality 

of Parallel Processing Elements" filed (DB00 1064-000, Micron no. 02-1269), and 

U.S. Patent Application Serial No. entitled "Method of Obtaining Interleave Interval 

for Two Data Values" filed (DB00 1065 -000, Micron no. 02-1268). 

BACKGROUND OF THE INVENTION 

[0002] The present invention relates generally to parallel processing and more particularly to 
balancing the work loads of the processing elements within a parallel processing system. 
[0003] Conventional central processing units ("CPU's"), such as those found in most 
personal computers, execute a single program (or instruction stream) and operate on a single 
stream of data. For example, the CPU fetches its program and data from a random access 
memory ("RAM"), manipulates the data in accordance with the program instructions, and 
writes the results back sequentially. There is a single stream of instructions and a single 
stream of data (note: a single operation may operate on more than one data item, as in X = Y 
+ Z, however, only a single stream of results is produced). Although the CPU may determine 
the sequence of instructions executed in the program itself, only one operation can be 
completed at a time. Because conventional CPUs execute a single program (or instruction 
stream) and operate on a single stream of data, conventional CPUs may be referred to as a 
single-instruction, single data CPU or an SISD CPU. 
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[0004] The speed of conventional CPUs has dramatically increased in recent years. 
Additionally, the use of cache memories enables conventional CPUs faster access to the 
desired instruction and data streams. However because conventional CPUs can complete only 
one operation at a time, conventional CPUs are not suitable for extremely demanding 
applications having large data sets (such as moving image processing, high quality speech 
recognition, and analytical modeling applications, among others). 

[0005] Improved performance over conventional SISD CPUs may be achieved by building 
systems which exhibit parallel processing capability. Typically, parallel processing systems 
use multiple processing units or processing elements to simultaneously perform one or more 
tasks on one or more data streams. For example in one class of parallel processing system, 
the results of an operation from a first CPU are passed to a second CPU for additional 
processing, and from the second CPU to another CPU, and so on. Such a system, commonly 
known as a "pipeline", is referred to as a multiple-instruction, single-data or MISD system 
because each CPU receives a different instruction stream while operating on a single data 
stream. Improved performance may also be obtained by using a system which contains many 
autonomous processors, each running its own program (even if the program running on the 
processors is the same code) and producing multiple data streams. Systems in this class are 
referred to as a multiple-instruction, multiple-data or MIMD system. 
[0006] Additionally, improved performance may be obtained using a system which has 
multiple identical processing units each performing the same operations at once on different 
data streams. The processing units may be under the control of a single sequencer running a 
single program. Systems in this class are referred to as a single-instruction, multiple data or 
SIMD system. When the number of processing units in this type of system is very large (e.g., 
hundreds or thousands), the system may be referred to as a massively parallel SIMD system. 
[0007] Nearly all computer systems now exhibit some aspect of one or more of these types 
of parallelism. For example, MMX extensions are SIMD; multiple processors (graphics 
processors, etc) are MIMD; pipelining (especially in graphics accelerators) is MISD. 
Furthermore, techniques such as out of order execution and multiple execution units have 
been used to introduce parallelism within conventional CPUs as well. 
[0008] Parallel processing is also used in active memory applications. An active memory 
refers to a memory device having a processing resource distributed throughout the memory 
structure. The processing resource is most often partitioned into many similar processing 
elements (PEs) and is typically a highly parallel computer system. By distributing the 
processing resource throughout the memory system, an active memory is able to exploit the 
very high data bandwidths available inside a memory system. Another advantage of active 
memory is that data can be processed "on-chip" without the need to transmit the data across a 
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system bus to the CPU or other system resource. Thus, the work load of the CPU may be 
reduced to operating system tasks, such as scheduling processes and allocating system 
resources. 

[0009] A typical active memory includes a number of interconnected PEs which are capable 
of simultaneously executing instructions sent from a central sequencer or control unit. The 
PEs may be connected in a variety of different arrangements depending on the design 
requirements for the active memory. For example, PEs may be arranged in hypercubes, 
butterfly networks, one-dimensional strings/loops, and two-dimensional meshes, among 
others. 

[0010] In typical active memories, load imbalances often occur such that some PEs are idle 
(i.e., without assigned tasks) while other PEs have multiple tasks assigned. To maximize the 
effectiveness of the active memory, it is desirable to balance the work load across all of the 
PEs. For example in an active memory having a multitude of identical PEs, it is desirable that 
each PE be assigned the same number of instructions by the central sequencer, thus 
maximizing the resources of the active memory. Additionally in an active memory having 
non-identical PEs, it may be desirable to assign more tasks to the PEs with greater processing 
capabilities. By balancing the load, the amount of time that one or more PEs is idle while* 
waiting for one or more other PEs to complete their assigned tasks is minimized. 
[0011] Thus, there exists a need for a method for balancing the load of a parallel processing 
system such that the resources of the parallel processing system are maximized. More 
specifically, there exists a need for a method for balancing the load of an active memory such 
that the resources of the active memory are maximized. 

SUMMARY OF THE INVENTION 

[0012] One aspect of the present invention relates to a method for balancing the load of a 
parallel processing system having a plurality of parallel processing elements (PEs) linked 
serially in a line with first and second ends, wherein each of the plurality of PEs has a local 
number of tasks associated therewith. The method comprises determining a total number of 
tasks present on the line, notifying each of the plurality of PEs of the total number of tasks, 
calculating a local mean number of tasks for each of the plurality of PEs, and calculating a 
local deviation for each of the plurality of PEs. The method also comprises determining a 
first local cumulative deviation for each of the plurality of PEs, determining a second local 
cumulative deviation for each of the plurality of PEs, and redistributing tasks among the 
plurality of PEs in response to the first local cumulative deviation and the second local 
cumulative deviation. 
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[0013] Another aspect of the present invention relates to method for assigning tasks 
associated with a parallel processing element (PE r ) serially connected to at least one other 
parallel processing element in a line of parallel processing elements, the PE r having a local 
number of tasks associated therewith. The method comprises notifying the PE r of a total 
number of tasks associated with the line, determining the PE r 's associated share of the total 
number of tasks, determining a first local cumulative deviation for the PE r , determining a 
second local cumulative deviation for the PE r , and redistributing tasks with at least one other 
connected parallel processing element in response to at least one of the first local cumulative 
deviation and the second local cumulative deviation. 

[0014] The present invention enables tasks to be distributed along a group of serially 
connected PEs so that each PE typically has X number of tasks or (X+l) number of tasks to 
perform in the next phase. The present invention may be performed using the hardware and 
software (i.e., the local processing capability) of each PE within the array. Those advantages 
and benefits, and others, will become apparent from description of the invention below. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] To enable the present invention to be easily understood and readily practiced, the: 
present invention will now be described for purposes of illustration and not limitation, in 
connection with the following figures wherein: 

[0016] FIG. 1 is a block diagram illustrating an active memory according to an embodiment 
of the present invention. 

[0017] FIG. 2 is a block diagram of a processing element for the active memory illustrated in 
FIG. 1 according to an embodiment of the present invention. 

[0018] FIG. 3 illustrates an array of processing elements as illustrated in FIG. 2 arranged in a 
line according to an embodiment of the present invention. 

[0019] FIG. 4 illustrates an operational process for balancing the load of the array of 
processing elements illustrated in FIG. 3 according to an embodiment of the present 
invention. 

[0020] FIGS. 5a - 5e are graphical representations of the load balancing method as applied 
to the array of processing elements illustrated in FIG. 3 according to an embodiment of the 
present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

[0021] As discussed above, parallel processing systems may be placed within one or more 
classifications (e.g., MISD, MIMD, SIMD, etc.). For simplicity, the present invention is 
discussed in the context of a SIMD parallel processing system. More specifically, the present 
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invention is discussed in the context of a SIMD active memory. It should be noted that such 
discussion is for clarity only and is not intended to the limit the scope of the present invention 
in any way. The present invention may be used for other types and classifications of parallel 
processing systems. 

[0022] FIG. 1 is a block diagram illustrating an active memory 10 according to an 
embodiment of the present invention. It should be noted that the active memory 10 is only 
one example of a device on which the methods of the present invention may be practiced and 
those of ordinary skill in the art will recognize that the block diagram of FIG. 1 is an 
overview of an active memory device 10 with a number of components known in the art being 
omitted for purposes of clarity. 

[0023] Active memory 10 is intended to be one component in a computer system. 
Processing within active memory 10 is initiated when the active memory 10 receives 
commands from a host processor (not shown), such as the computer system's CPU. A 
complete processing operation (i.e., data movement and processing) in the active memory 10 
may consist of a sequence of many commands from the host to the active memory device 10. 
[0024] Active memory 10 is comprised of a host memory interface ("HMI") 12, a bus 
interface 14, a clock generator 16, a task dispatch unit ("TDU") 18, a DRAM control unit 
("DCU") 20, a DRAM module 22, a programmable SRAM 24, an array control sequencer 26, 
and a processing element array 28, among others. 

[0025] The HMI 12 provides an input/output channel between the host (such as a CPU, not 
shown) and the DRAM module 22. In the current embodiment, the HMI 12 receives 
command (cmd), address (addr), and data signals (among others) from and sends data and 
ready (rdy) signals (among others) to the host. The HMI 12 approximates the operation of a 
standard non-active memory so that the host, without modifications, is compatible with the 
active memory 10. 

[0026] The HMI 12 may be similar in its operation to the interface of a synchronous DRAM 
as is know in the art. Accordingly, the host must first activate a page of data to access data 
within a DRAM module 22. In the current embodiment, each page may contain 1024 bytes of 
data and there may be 16,384 pages in all. Once a page has been activated, it can be written 
and read through the HMI 12. The data in the DRAM module 22 may be updated when the 
page is deactivated. The HMI 12 also sends control signals (among others) to the DCU 20 
and to the processing element array 28 via the task dispatch unit 18. 

[0027] The HMI 12 may operate at a frequency different than that of the frequency of the 
master clock. For example, a 2x internal clock signal from clock generator 16 may be used. 
Unlike a traditional DRAM, the access time for the HMI 12 uses a variable number of cycles 
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to complete an internal operation, such as an activate or deactivate. Thus the ready signal 
(rdy) is provided to allow the host to detect when a specific command has been completed. 
[0028] The bus interface 14 provides and input/output channel between the host and the 
TDU 18. For example, the bus interface 14 receives column select (cs), write command (w), 
read command (r), address (addr), and data signals (among others) from and places interrupt 
(intr), flag, and data signals (among others) onto the system bus (not shown). The bus 
interface 14 also receives signals from and sends signals to TDU 18. 

[0029] The clock generator 16 is operable to receive an external master clock signal (xl) and 
operable to provide the master clock signal (xl) and one or more internal clock signals (x2, 
x4, x8) to the components of the active memory. It should be apparent to one skilled in the 
art that other internal clock signals may be produced by the clock generator 16. 
[0030] The TDU 18 communicates with the bus interface 14, the HMI 12, the programmable 
SRAM 24, the array control sequencer 26, and the DCU 20. In the current embodiment, the 
TDU 1 8 functions as an interface to allow the host to issue a sequence of commands to the 
array control sequencer 26 and the DCU 20. Task commands from the host may be buffered 
in the TDU's FIFO buffers to allow a burst command to be issued. Commands may contain 
information on how the tasks in the array control sequencer 26 and the DCU 20 should be 
synchronized with one another, among others. 

[0031] The DCU 20 arbitrates between the TDU 1 8 and the HMI 12 and sends commands to 
the DRAM modules 22 and the processing element array 28. The DCU 20 also schedules 
refreshes within the DRAM modules 22. In one embodiment, the DRAM modules 22 of the 
active memory 10 may be comprised of sixteen 64k xl28 eDRAM (or embedded DRAM) 
cores. Each eDRAM core may be connected to an array of sixteen PEs, thus providing 256 
(16 x 16) PEs in all. 

[0032] The programmable SRAM 24 functions as a program memory by storing commands 
issued by the TDU 18. For example, the TDU 18 may transmit a "write program memory 
acldress" command which sets up a start address for a write operation and a "write program 
memory data" command which writes/a memory location and increments the program 
memory write address, among others. The programmable SRAM 24, in the current 
embodiment, has both an address register and a data output register. 
[0033] The array control sequencer 26 may be comprised of a simple 16 bit minimal 
instruction set computer (16-MISC). The array control sequencer 26 communicates with the 
TDU 18, the programmable SRAM 24, and the DCU 20, and is operable to generate register 
file addresses for the processing element array 28 and operable to sequence the array 
commands, among others. 
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[0034] The processing element array 28 is comprised of a multitude of processing elements 
("PEs") 30 (see FIG. 2) connected in a variety of different arrangements depending on the 
design requirements for the processing system. For example, processing units may be 
arranged in hypercubes, butterfly networks, one-dimensional lines/loops, and two- 
dimensional meshes, among others. In the current embodiment (e.g., as shown in FIG. 3), the 
processing elements 30 are arranged in a line (i.e., a one-dimensional array). The processing 
element array 28 communicates with the DRAM module 22 and executes commands received 
from the programmable SRAM 24, the array control sequencer 26, the DCU 20, and the HMI 
12. Each PE in the processing element array 28 includes dedicated H-registers for 
communication with the HMI 12. Control of the H-registers is shared by the HMI 12 and the 
DCU 20. It should be noted that "line" refers to at least two serially connected PEs. It should 
be apparent to one skilled in the art that "line" is intended to include PE's arranged in a linear 
array (e.g., rows, columns, diagonals, etc.) and other non-linear shapes. It should also be 
apparent to one skilled in the art that serially connected arrays having uniform and/or varied 
distances between one or more of the PEs are within the scope of the present invention. 
[0035] Referring now to FIG. 2, a block diagram of a PE 30 according to one embodiment of 
the present invention is illustrated. PE 30 includes an arithmetic logic unit ("ALU") 32, Q- 
registers 34, M-registers 36, a shift control and condition register 38 (also called "condition 
logic" 38), a result register pipeline 40, and register file 42. The PE 30 may also contain other 
components such as multiplexers 48 and logic gates (not shown), among others. 
[0036] In the current embodiment, the Q-registers 34 are operable to merge data into a 
floating point format and the M-Registers 36 are operable to de-merge data from a floating, 
point format into a single magnitude plus an exponent format. The ALU 32 is a multiplier-, 
adder operable (among others) to receive information from the Q-registers 34 and M-registers 
36, execute tasks assigned by the TDU 18 (see FIG. 1), and transmit results to the shift 
control and condition logic 38 and to the result register pipeline 40. The result register 
pipeline 40 is operable to communicate with the register file 42, which holds data for transfer 
into or out of the DRAM modules 22 via a DRAM interface 44. Data is transferred between 
the PE and the DRAM module 22 via a pair a registers, one register being responsive to the 
DCU 20 and the other register being responsive to the PE 30. The DRAM interface receives 
command information from the DCU 20. The DRAM interface 44 also permits the PE 30 to 
communicate with the host through the host memory access port 46. 

[0037] In the current embodiment, the H-registers 42 are comprised of synchronous SRAM 
and each processing element within the processing element array 28 contains eight H-registers 
42 so that two pages can be stored from different DRAM locations, thus allowing the 
interleaving of short i/o bursts to be more efficient. Result register pipeline 40 is also 

7 

0O435057.DOC 



connected to one or more neighborhood connection registers ("X-register") (not shown). The 
X-register links one PE 30 to its neighboring PE's 30 in the processing element array 28. 
[0038] The reader desiring more information about the hardware shown in FIGs. 1 and 2 is 
directed to UK Patent application (serial no. not yet assigned) entitled "Control of Processing 
Elements in Parallel Processors" filed 17 September 2002, (Micron no. 02-1604) which is 
hereby incorporated by reference. Details about the PEs may also be found in UK Patent 
Application No. 021562.2 entitled "Host Memory Interface for a Parallel Processor" filed 17 
September 2002, (Micron no. 02-0703) which is hereby incorporated by reference. 
[0039] FIG. 3 is a simplified diagram showing the interconnections of an array of the 
processing elements 30 (as illustrated in FIG. 2) arranged in a line 50 according to an 
embodiment of the present invention. In the current embodiment, line 50 is comprised of 
eight (8) PEs 30 (i.e., PEo, PE,, . . . PE 7 ) which are interconnected via their associated X- 
register links. It should be noted that the number of PEs 30 included in line 50 may be altered 
while remaining within the scope of the present invention. As illustrated in FIG. 3, PEj 
through PE 6 are operable to communicate with both their left and right neighbors. For 
example, PE! can communicate with PE 0 (i.e., PE/s left neighbor) and with PE2 (i.e., PE^s 
right neighbor). In the current embodiment, the line's 50 left end PE (i.e., PE 0 ) is operable to 
communicate with it's right neighbor (i.e., PEi). Similarly, the line's 50 right end PE (i.e., 
PE 7 ) is operable to communicate with it's left neighbor (i.e., PE 6 ). It should be noted, 
however, that each end PE (i.e., PE 0 and PE 7 ) may also be operable to communicate with a PE 
from another line (e.g., to link two or more lines in an array) or to communicate with each 
other (e.g., to permit a wrap function). In the current embodiment, every PE 30 in the line 50 
receives instructions from a single TDU 18 as discussed in conjunction with FIG. 1. 
[0040] FIG. 4 illustrates an operational process 60 for balancing the work loads between the 
PEs 30 in line 50 according to an embodiment of the present invention. Operational process 
60 begins by determining the total number of tasks (V) present for the line in operation 61. 
Each PE r (where r = 0 to 7, e.g., PE 0 , PE 1? . . . PE 7 ) in the line 50 has a local number of tasks 
(v r ) associated therewith. 

[0041] In the current embodiment, a partial sum (S r ) of the number of tasks (v r ) is passed 
from a first end to a second end (e.g., from the leftmost PE to the rightmost PE). The sum 
(S r ) is initialized to the number of tasks on the first end (here 5 0 = v 0 ). The partial sum (S r ) is 
then passed serially from PEr to PE^i (i.e., from left to right). Each PEr receives the partial 
sum and adds it's own number of tasks to the partial sum (i.e., S r =S r _/ + v r ). After adding it's 
number of tasks (v r ), PEr passes the partial sum (S r ) to PE,^. When the partial sum (S r ) 
reaches the right end (i.e., PE N _i), the last tasks (v M/ ) are added to form the total sum of values 

8 

0O435O57.DOC 



(V). The sum (V) can be expressed by the equation V = v i ;j where //represents the 

number of PEs 30 in the line 50, and v f - represents the local number of tasks associated with a 
local PE r in the line 50. It should also be noted that in the current discussion, "local" refers to 
the values or functions associated with a single PE within the line, whereas "global" refers to 
the values or functions associated with the entire line of PEs. 

[0042] Referring to FIG. 3, PEo has 3 tasks associated therewith (i.e., v 0 =3). Thus, the initial 
value for the partial sum is also 3 (i.e., S 0 = v 0 =3). The partial sum is then passed to PEi, 
which has 6 tasks associated therewith (i.e., v ; =6). The number of tasks associated with PE r 
is then added to the partial sum S 0 (i.e., 5/ = S 0 + vi) to obtain S f = 9. The partial sum 
continues to be passed serially until it reaches the right end of line 50 (i.e., PE 7 ). As 
illustrated in FIG. 3, the total sum of tasks on line 50 is equal to 40 (i.e., V— 40). 
[0043] After the total number of tasks (V) present on the line is determined in operation 61 , 
each PE 30 is notified of the total number of tasks (V) on the line 50 by operation 62. In the 
current embodiment, the total number of tasks (V) is transferred from right to left (i.e., from 
PE 7 to PE 0 ). Each PE in line 50 locally saves the total number of tasks (V). It should be 
apparent to one skilled in the art that the direction of travel in operations 61 and 62 may be 
altered while remaining within the scope of the present invention. For example, the direction 
of travel for operations 61 and 62 may be changed from left/right to right/left and from 
right/left to left/right, respectively. 

[0044] After the total number of tasks (V) present on the line 50 is distributed in operation 
62, the local mean number (M r ) of tasks for each PE r is calculated in operation 63. In the 
current embodiment, the local mean value is computed using the rounding function 
M r = Trunc((V + E r )l N) (where M r represents the local mean for PE r , iV represents the 
total number of PEs 30 in the line 50, and E r represents a number in the range of 0 to (Af-1)) 
to ensure that no instructions are lost or "gained" during the rounding process if the value of V 

+ Nis not an integer (i.e., to ensure that V = ^ A M t . , where N represents the number of PEs 

i=0 

30 in the line 50, and M, represents the local mean of tasks associated with a local PE r in the 
line 50). The rounding function is discussed in more detail in U.S. Patent Application Serial 

No. entitled "Rounding Algorithm for a Plurality of Parallel Processing Elements" 

filed (DB00 1064-000, Micron no. 02-1269) and incorporated in its entirety by 

reference herein. 

[0045] For example assuming that 43 tasks (i.e., V= 43) were to be shared by eight PEs in a 
line (i.e., PEo through PE 7 ), and without using the rounding function, the local mean for each 
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PE would be PE r = 5.375 (i.e., 43 -s- 8 = 5.375). If the fraction 43/8 is designated to round 
down for each PE (i.e., 43 8 = 5), then the sum of the means for all of the individual PEs 
(i.e., PEo through PE 7 ) will equal to 40. Thus, three (43 - 40 = 3) tasks are lost. In contrast, if 
the fraction 43/8 is designated to round up for each PE (i.e., 43 8 = 6), then the sum of the 
means for all of the individual PEs (i.e., PEo through PE 7 ) will equal to 48. Thus, five (48 - 
43 = 5) tasks are gained. 

[0046] The rounding function M r = Trunc((V + E r )l N) prevents tasks from being lost or 
gained. In the current embodiment, each PE is assigned a different E r value for controlling 
the rounding. The simplest form for the function E is the case in which E r = P r the number of 
the PE. For example, for PE 0 , E 0 = 0; for PEi, Ej = 1; for PE 2 , E 2 = 2; etc. By assigning each 
PE 30 a different E r value, the rounding function can be controlled such that some of the local 
means are rounded up and some of the local means are rounded down, thus insuring that 

V = y^M, . It should be noted that in the current embodiment, the local mean for each PE 

30 in line 50 is computed in parallel with the local means of the other PEs in the line 50. 
[0047] Table 1 , for example, illustrates the local mean calculation for the case in which the 
total number of tasks on the line 50 is equal to forty-three (43). Referring to Table 1, it is 
apparent that the rounding function controls the rounding such that M 0 through M 4 are 
rounded to five (5), whereas M s through M 7 are rounded to six (6). The sum of the values of 
M 0 through M 7 is equal to 43, which is equal to the total number of tasks (V) on the line 50. 
Thus, tasks were neither lost nor gained due to rounding. 
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Table #1 - Local Mean Calculation For V= 43, N= 8. 



[0048] Table 2 illustrates the local mean calculation for the line of PEs as shown in FIG. 3. 
Each PE r calculates its local mean (M r ). As previously stated in conjunction with operation 
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61 , the total number of tasks (V) associated with line 50 is equal to 40. Thus, each PE r will 
calculate a local mean equal to five (5). As expected, the sum of the local means (M r ) is equal 
to the total number of tasks (V). 
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Table #2 - Local Mean Calculation For K= 40, N= 8. 



[0049] After the local means (M r ) are computed in operation 63, the local deviation D r is 
calculated for each PE r in operation 64. In the current embodiment, the local deviation is 
simply the difference between the local value and the local mean (i.e., D r = v r - Mi). 
[0050] After the local deviations (Z) r ) are computed in operation 64, a first local cumulative 
deviation (L r ) for each PE is determined in operation 65. In the current embodiment, the first 
local cumulative deviation (L r ) is calculated from left to right (i.e., from PE 0 to PE 7 ). The first 
local cumulative deviation (L r ) is determined in a similar manner as that used to form the 
partial value sum (V) in operation 61. In operation 65, however, the local deviations (D r ) are 
summed instead of the number of tasks (v r ). 

[0051] Additionally, a PE r does not need to be notified of the local deviations for 
"downstream" PEs (i.e., PEs which come after PE r depending on the direction of travel). In 
contrast, the first local cumulative deviation for a specific PE r is the sum of the local 
deviations for each of the upstream (i.e., previous) PEs within the line. For example in the 
current embodiment in which the direction of travel is left to right, the first local cumulative 
deviation for PEj is L 3 = D 0 + £>/ + D 2 . Thus, PEj does not need to be notified of the first 
local cumulative deviations for the downstream PEs (i.e., for PE, through PE 7 ). The first 

i=r-J 

local cumulative deviation (L r ) can be represented by the equation: L r = ^Z) f , where r 
represents the local number of the specific PE for which the first local cumulative deviation is 
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being calculated, and A represents the local deviation associated with the "upstream" PEs in 
the line 50. 

[0052] After first local cumulative deviation (L r ) is determined in operation 65, the second 
local cumulative deviation (R r ) is determined for each PE in operation 66. In the current 
embodiment, the second local cumulative deviation (R T ) is simply evaluated by observing that 
the sum of deviation across the whole line 50 must be zero. Thus, the second local 
cumulative deviation can be expressed as R r = ~D r - L r . In the current embodiment, R r = -D r - 
L r is evaluated locally, in parallel on all PEs in the line. It should be noted, however, that R r ~ 
-D r - L r may also be evaluated serially while remaining within the scope of the present 
invention. 

[0053] After the first local cumulative deviation (Z r ) and the second local cumulative 
deviation (R r ) are determined for each PE in operations 65 and 66, respectively, the tasks 
associated with each PE are redistributed among the PEs relative to the values obtained for L r 
and R r . In the current embodiment, a positive L r parameter represents the number of tasks, to 
the left of the current PE r , that are "surplus to the requirements". A negative L r value, 
represents the "shortfall" in the number of tasks to the left to the current PE r . Thus if L r is 
positive, this number of tasks will be received from the left PE, whereas if L r is negative, this 
number of tasks will need to be transmitted to the left PE. Similarly, if R r is positive, this 
number of tasks will be received from the right PE, whereas if R r is negative, this number of 
tasks will need to be transmitted to the right PE. 

[0054] If the local deviation (JD r ) is negative, one or more of the received values will be 
"absorbed" by the local PE to make up the local deficit. The other received values (if any)- 
will be transmitted, either from left to right, or from right to left depending on the values for 
L r and R r . On occasion, some PEs may start off with no values at all, these PEs may have to 
"mark time" until they receive a value. It should be noted that in the current embodiment, 
only one task (per direction) is transferred between PEs at a time. For example, only one task 
may be transferred between PE 3 and PE 2 , however, one task may be transferred between PE 3 
and PE 4 at the same time. After each successful transmission or receipt, the values for D FJ L n 
and R r are updated. The redistribution stage only terminates when D r = L r = R r = 0 for all 
PEs. 

[0055] FIGS. 5a - 5e are a graphical representation of the load balancing method as applied 
to the line 50 of processing elements illustrated in FIG. 3 according to an embodiment of the 
present invention. In FIGS. 5a - 5e, each PE (i.e., PE 0 . . . PE 7 ) is represented by a bar, the 
height of which indicates the number of tasks (v r ) associated with the specific PE. FIGS. 5a - 
5e also include the local deviation (D r ), the first local cumulative deviation (L r ), and the 
second local cumulative deviation (R r ) for each PE r within the line 50. Referring to FIG. 5a 
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for example, PE 3 has a local deviation of +2 (i.e., D 3 = 2), a first local cumulative of deviation 
- 4 (i.e., L 3 = -4), and a second local cumulative deviation +2 (i.e., R 3 - +2). 
[0056] FIG. 5a represents the line 50 as illustrated in FIG. 3 (i.e., the initial number of tasks 
assigned to each PE in FIG. 5 is the same as illustrated in FIG. 3). Accordingly, the total 
number of tasks (V) is calculated to be equal to forty (40). As discussed in conjunction with 
Table 2, the local mean number (M r ) is equal to five (5) for each PE. In the current 
embodiment, the local deviation (D r ), the first local cumulative deviation (L r ), and the second 
local cumulative deviation (R r ) are calculated as discussed in conjunction with FIG. 4. 
[0057] FIG. 5b represents the line 50 after a first redistribution of tasks has taken place. It 
should be noted that the number of tasks transferred and the direction in which the tasks are 
transferred are related to the first local cumulative deviation (Z r ) and the second local 
cumulative deviation (R r ). As discussed above in conjunction with FIG. 4, a positive L r 
represents the number of tasks that will be received by the local PE r from its left-neighboring 
PE r ./, whereas a negative L r represents the number of tasks that the local PE r transmits to its 
left-neighboring PE r _/. Similarly, a positive R r represents the number of tasks that will be 
received by the local PE r from its right-neighboring PE r+/ , whereas a negative R r represents 
the number of tasks that the local PE r transmits to its right-neighboring PE r+/ . It should 
further be noted that in the current embodiment, only one task is transferred between PEs at a 
time. 

[0058] PE 3 , for example, initially has a first local cumulative deviation of - 4 (i.e., L 3 = -4), 
and a second local cumulative deviation of +2 (i.e., R 3 = +2). Accordingly, one task is 
transferred to its left-neighboring PE r . y (i.e., PE 2 ) and one task is received from its right- 
neighboring PE r+/ (i.e., PE 4 ) (as indicated by the arrows in FIG. 5b). After the transfer is 
completed, the local deviation, the first local cumulative deviation and the second local 
cumulative deviation are recalculated for each PE, for example, D 3 = +2, L 3 = -3 and R 3 = +1, 
respectively, after the first transfer. 

[0059] FIG. 5c illustrates the next transfer of tasks between the PEs. PE 3 , for example, 
transfers another task to its left-neighboring PE r . y (i.e., PE 2 ) and receives another task from its 
right-neighboring PE r+/ (i.e., PE 4 ). After the transfer is completed, the local deviation, the 
first local cumulative deviation and the second local cumulative deviation are recalculated for 
each PE, for example, D 3 = +2, L 3 = -2 and R 3 = 0, respectively. 

[0060] FIG. 5d illustrates the next transfer of tasks between the PEs. PE 3 , for example, 
transfers another task to its left-neighboring PE r ./ (i.e., PE 2 ). However because R 3 -0 (and 
likewise L 4 = 0), PE 3 does not receive another task from its right-neighboring PE r+/ (i.e., PE 4 ). 
After the transfer is completed, the local deviation, the first local cumulative deviation and the 
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second local cumulative deviation are recalculated for each PE, for example, D 3 = +1, L 3 = -1 
and R 3 = 0, respectively. 

[0061] FIG. 5e illustrates the final transfer of tasks between the PEs. As seen in FIG. 5d, 
only PE 3 and PE 2 have non-zero local and cumulative deviation values. Thus, tasks will be 
transferred only between PE 3 and PE 2 . Referring to FIG. 5e, PE 3 transfers another task to its 
left-neighboring PE r _ 7 (i.e., PE 2 ). Again because R 3 = 0 (and likewise L 4 = 0), PE 3 does not 
receive another task from its right-neighboring PE r+/ (i.e., PE 4 ). After the transfer is 
completed, the local deviation, the first local cumulative deviation and the second local 
cumulative deviation are recalculated for each PE. As illustrated, the local and cumulative 
deviations for all PEs equals zero. Thus, the load is balanced among all PEs and further 
transfers are not necessary. 

[0062] It should be recognized that the above-described embodiments of the invention are 
intended to be illustrative only. Numerous alternative embodiments may be devised by those 
skilled in the art without departing from the scope of the following claims. For example, the 
present invention may form an elemental part of a larger process such as for a set of 
processors that are connected in a line, but may also be part of a rectangular array and cuboid 
array, among others. 
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